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Abstract 

Linear carbon-chain radicals C„H for n = 2-7 have been studied with correlation consistent valence and core-valence 
basis sets and the coupled cluster method RCCSCKT). Equilibrium structures, rotational constants, and dipole moments are 
reported and compared with available experimental data. The ground state of the even-* series changes from to “FI as 
the chain is extended. For C 4 H, the 2 X + state was found to lie only 72 cm -1 below the 2 n state in the estimated complete 
basis set limit for valence correlation. The C 2 H and C 3 H anions have also been characterized. 


1. Introduction 

The acetylenic linear carbon-chain radicals C„H 
for n — 2 to 6 have all been detected in the interstel- 
lar medium (see citations in Refs. [1-3]), with C 2 H 
particularly plentiful. As reactive precursors in 
chain-lengthening astrochemical processes, they and 
the related cyanopolyyne C n N radicals play impor- 
tant roles in the evolution of molecular species within 
dense interstellar and circumstellar molecular clouds. 
The members of the odd-n series all have 2 U ground 
states with characteristically high dipole moments, 
which make them exceptional candidates for astro- 
nomical detection. For the even -n series, there are 
low-lying states of 2 S + and 2 II symmetries, where 
the 2 S, + state has a much smaller dipole moment 
than the 2 fl state, though large enough to be de- 


1 E-mail: woon@heda.molres.org 


tected. The ordering of the states unquestionably 
changes as the length of the chain grows. There is no 
dispute that the ground state of C 2 H (ethynyl) is a 
2 2 + , while that of C 6 H (hexatriynyl) is 2 IT. How- 
ever, the ground state of the C 4 H butadiynyl radical 
lies near the point where the ordering of the two 
states reverses. The experimental spectra of C 4 H are 
consistent with a 2 S + ground state [4-7], but the 
state separation has not been measured. Calculations 
by Horn and Botschwina included in Ref. [7] found 
the 2 £ + state to be more stable, while Natterer and 
Koch [8] asserted that the ground state is of ~ II 
symmetry. Their work, however, left some points 
unsettled. Thus one goal of the present study is to 
reexamine this issue in order to ascertain if the 
proper ground state assignment has been made for 
C 4 H, as well as to compute accurate values for the 
state separations of C 2 H and C 6 H. 

Although each of the C„H species containing up 
to ten carbon atoms has been the focus of previous 
ab initio studies, there have been only three efforts to 
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characterize all or part of the sequence at the same 
level of theoretical treatment. Cooper and Murphy 
[1] reported the results of RHF calculations on the 
odd-rt C 3 H, C 5 H, and C 7 H species using very lim- 
ited basis sets. Subsequently, Pauzat et al. [3] studied 
both even-and odd-/i radicals from C 2 H through 
C 10 H at the UHF/3-21G and UHF/SVP levels. 
Finally, the work by Natterer and Koch [8] men- 
tioned above considered the even -n radicals C 2 H, 
C 4 H and C 6 H. Only the last paper employed corre- 
lated methods. The complete series has not been 
systematically investigated with correlated methods 
and robust basis sets. Therefore, the second goal of 
this work is to study C 2 H through C 7 H with state- 
of-the-art ab initio theory and substantial basis sets 
in order to predict molecular properties including 
structures, rotational constants, and dipole moments. 

Finally, the electron affinities (EA) of C 2 FI and 
C 3 H, the two C„H species for which experimental 
values have been measured, were also investigated in 
this work. The C 2 H _ and C 3 H” anions are of 
potential interest in astrochemistry due to their stabil- 
ity and the well-known importance of ion-molecule 
reactions under astrophysical conditions. Neither EA 
has received much theoretical attention. In particular, 
C 3 H“ has only been studied in one previous effort 
[9], which utilized a moderately-sized basis set (6-31 
+ + G * * ) and limited recovery of correlation (CID) 
but did not report the EA. It did, however, provide 
the first evidence that the geometry of the ground 
state of C 3 H~ is not linear and that there is also an 
intriguing bound excited state. By using larger basis 
sets and a correlation method [RCCSD(T)] that per- 
forms well at describing EAs, we can improve upon 
earlier predictions. 


2. Methodology 

The correlation consistent basis sets (cc-pVXZ) of 
Dunning and co-workers [10-12] have been utilized 
in this work. An extensive series of benchmark 
calculations [13-18] on small systems have amply 
demonstrated the well behaved convergence associ- 
ated with the use of these basis sets, which makes it 
possible to estimate complete basis set (CBS) limits 
for various properties when sets of at least double 
through quadruple zeta quality can be used. Even in 


cases where the larger sets cannot be used, there are 
systematic trends that may be exploited in order to 
improve the accuracy of the ab initio predictions. 

The primary correlation method employed in this 
study was the restricted coupled-cluster method 
RCCSEKT) [19] as implemented in the MOLPRO 
suite of codes 2 . An RHF zeroth-order wave func- 
tion was used. For the species with 2 n states, the 
two degenerate components were state-averaged to 
provide orbitals with the proper symmetry. In order 
to evaluate the performance of RCCSEKT) for de- 
scribing the 2 S^- 2 I1 state separations of the even-« 
series, internally contracted multireference configura- 
tion interaction (CMRCI) calculations [20] were car- 
ried out for C 2 H using a complete active space 
self-consistent field (CASSCF) reference function 
[21]. The multireference Davidson correction 
(CMRCI + 0) was also applied [22]. Although most 
of the correlation calculations were restricted to the 
valence electrons, the structures of C 2 H and C 3 H 
were reoptimized with the new cc-pCVTZ set for 
carbon [12], which includes additional optimized 
functions located in the core and core-valence re- 
gion. 

Dipole moments at the RCCSEKT) level were 
determined using the finite field approximation. The 
basis sets were enlarged to include extra diffuse 
functions by using the aug-cc-pVXZ sets [11] of 
double or triple zeta quality. 


3. Structures, rotational constants, and dipole mo- 
ments of C n H species 

Optimized equilibrium bond lengths and the asso- 
ciated rotational constants for both even-and odd-n 
C n H radicals determined at the RCCSEKT) level of 
theory are reported in tblrl. The largest valence sets 
used for each species were dictated by available 
computational resources. The cc-pVQZ sets were 
used for the full geometry optimizations of C 2 H and 
C 3 H and for single-point calculations on C 4 H. The 


2 MOLPRO is a package of ab initio programs written by H.-J. 
Werner and P.J. Knowles with contributions by J. AJmlof, R.D. 
Amos, M.J.O. Deegan, S.T. Elbert, C. Hampel, W. Meyer, K.A. 
Peterson, R.M. Pitzer, A.J. Stone and P.R. Taylor. 
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cc-pVTZ sets could be used for C 5 H and single 
points for C 6 H, while only the cc-pVDZ sets could 
be used for optimizing C 7 H. 

One benefit in examining a large number of 
species simultaneously is that trends may emerge 
that are not apparent when examining individual 
molecules. Here there are a total of 9 C-H and 30 
C-C bonds. The latter vary in character from single 
to triple bonds. When the basis set is improved from 
cc-pVDZ to cc-pVTZ quality, both C-H and C-C 
bond lengths decrease, but the effect is remarkably 
uniform not just for very different C-C bonds, but 
also for C-C versus C-H bonds. The average ratio 
( r e (cc-pVTZ) /r e (cc-pVDZ)) for 6 C-H bonds is 
0.9855, but the high-low range is only 0.9853 to 
0.9856. The range of 0.9838 to 0.9871 for 13 C-C 
bonds is larger, but the average is again 0.9855. The 
predictive power of this parameter for estimating 
bond lengths at the cc-pVTZ level may be gauged by 


considering an example. For C 3 H ( 2 n), the average 
parameter of 0.9855 produces bond lengths that yield 
B t = 11.030 GHz, which is essentially the same as 
the optimized cc-pVTZ result of 11.028 GHz. If the 
low and high extremes are used, they lead to respec- 
tive B e values of 11.065 and 10.995 GHz, an inter- 
val of only 70 MHz. The average scaling parameter 
was use to estimate bond lengths for single point 
cc-pVTZ calculations for the two states of C 6 H in 
order to produce a more accurate prediction for the 
state separation (see Section 4). 

For the improvement from cc-pVTZ to cc-pVQZ 
quality, there are only 3 C— H and 4 C— C bond 
lengths upon which to base scaling factors. Here the 
average C-H and C-C scaling factors (r e (cc- 
pVQZ)/r e (cc-pVTZ)) are somewhat different, 0.9998 
for C-H bonds and 0.9970 for C-C bonds. The 
variation for C-H bonds is again very small, and the 
range for C-C bonds is only 0.9967 to 0.9978. The 


Table 1 

Bond lengths r c and rotational constants B c for the "'X* and 11 states of C„H (fl — 2—7) at the RCCSEKT) level of theory 
Species State Basis set r e (A) &e (GHz) £ 0 (GHz) 

H-C^i C,-C 2 C 2 -C 3 C 3 -C 4 C 4 - C, C 5 -Q C 6 -C 7 calc. exp. 


c 2 h 

2 X + 

cc-pVDZ 

1.0801 



cc-pVTZ 

1.0645 



cc-pVQZ 

1.0644 



est. CBS limit 

1.0644 


2 11 

cc-pVDZ 

1.0876 



cc-pVTZ 

1.0716 



cc-pVQZ 

1.0712 



est. CBS limit 

1.0712 

c 3 h 

2 n 

cc-pVDZ 

1.0813 



cc-pVTZ 

1.0654 



cc-pVQZ 

1.0653 



est. CBS limit 

1.0653 

c 4 h 

2 X + 

cc-pVDZ 

1.0793 



cc-pVTZ 

1.0635 



cc-pVQZ a 

1.0633 


2 n 

cc-pVDZ 

1.0805 



cc-pVTZ 

1.0648 



cc-pVQZ a 

1.0645 

c 5 h 

2 n 

cc-pVDZ 

1.0805 



cc-pVTZ 

1.0648 

c 6 h 

2 X + 

cc-pVDZ 

1.0796 



cc-pVTZ a 

1.0639 


2 n 

cc-pVDZ 

1.0804 



cc-pVTZ a 

1.0648 

c 7 h 

2 n 

cc-pVDZ 

1.0804 


1.2348 

1.2148 

1.2108 

1.2088 

1.3154 

1.2953 

1.2910 

1.2888 


1.2661 

1.3670 



1.2484 

1.3468 



1.2456 

1.3426 



1.2443 

1.3405 



1.2334 

1.3947 

1.2396 


1.2146 

1.3787 

1.2200 


1.2104 

1.3744 

1.2157 


1.2470 

1.3621 

1.3179 


1.2286 

1.3444 

1.2984 


1.2244 

1.3402 

1.2942 


1.2457 

1.3567 

1.2899 

1.3371 

1.2276 

1.3392 

1.2727 

1.3176 

1.2349 

1.3881 

1.2397 

1.3880 

1.2171 

1.3680 

1.2217 

1.3679 

1.2389 

1.3734 

1.2598 

1.3462 

1.2210 

1.3535 

1.2415 

1.3267 

1.2395 

1.3724 

1.2608 

1.3357 


42.341 

43.723 

43.972 

44.096 43.675 

37.917 
39.101 
39.329 
39.443 
10.712 
11.028 
11.085 

11.112 11.186 

4.578 
4.708 

4.738 4.759 

4.484 

4.612 

4.642 

2.302 

2.366 2.395 

1.2411 1.339 

1.2231 1.379 

1.3185 1.332 

1.2994 1.371 1.386 

1.2997 1.3288 0.842 


a Single-point calculation at geometry derived from smaller set(s) (see text). 
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two average scaling factors were used to estimate 
bond lengths in order to run the cc-pVQZ calcula- 
tions for C 4 H that appear in Table 1 and play an 
important role in Section 4. 

The effect of including core and core-valence 
correlation can be approximated by running all-elec- 
tron calculations at the cc-pCVTZ level and observ- 
ing the changes in bond lengths with respect to 
valence-electron calculations at the cc-pVTZ level. 
The structures of C 2 H (both states) and C 3 H were 
therefore reoptimized. The C-H bonds are essen- 
tially unchanged, but the C-C bonds decrease by an 
average of 0.0034 A, which is consistent with shifts 
observed in C-C bonds in C 2 [23], C 3 [24] and 
MgCCH [25]. 

Rotational constants are the principle point of 
comparison between the present calculations and ex- 
perimental or astronomical data. Measured B 0 val- 
ues are indicated in Table 1, as collected by Pauzat 
et al. [3]. Passing over C 2 H for the moment, the 
agreement is very good between the RCCSEKT) 
values (valence correlation, CBS limits or largest 
basis set result) and the empirical data for species 
where there is no dispute over the ground state 
assignment. The computed B t values for C 3 H, C 5 H 
and C 6 H fall only 74, 29 and 57 MHz below the 
measurements, where some adjustment for zero-point 
vibrational motion and other theoretical approxima- 
tions should be expected. The estimated cc-pVQZ 
geometries for C 4 H yield B c values that fall 21 and 
117 MHz below the result of Gottlieb et al. [5] for 
the and 2 I1 states, respectively, which offers 

tentative support for a 2 X + ground state. When 
core-valence contractions of 0.0034 A per C-C 
bond are added to the valence CBS limits or best 
basis set results noted above, the values of B e be- 
come 11.165, 4.763, 2.379 and 1.384 GHz for the 
ground states of C 3 H to C 6 H, respectively. The 
respective errors of 21, 4, 16 and 2 MHz are signifi- 
cantly reduced from the valence-only results. 

Theory overestimates B t for C 2 H by more than 
400 MHz. The RCCSD(T) results were tested against 
the CMRCI method. The CMRCI estimated CBS 
limits for the 2 S + and 2 II states were found to be 
44.003 and 39.372 GHz, respectively, which differ 
from RCCSD(T) by 93 and 71 MHz. The behavior 
of the RCCSD(T) method is not anomalous. The 
large difference between theory and experiment may 


Table 2 

Dipole moments ja c of the "X* and 2 n states of C„H at the 
RCCSIXT) level 


Species 

Basis set 


Mc(D) 


C 

H 

2 2 + 

2 n 

C 2 H 

cc-pVDZ 

cc-pVDZ 

0.766 

2.829 


aug-cc-pVDZ 

cc-pVDZ 

0.772 

2.965 


aug-cc-pVDZ 

aug-cc-pVDZ 

0.774 

2.968 


cc-pVTZ 

cc-pVTZ 

0.774 

2.960 


aug-cc-pVTZ 

cc-pVTZ 

0.771 

3.006 


aug-cc-pVTZ 

aug-cc-pVTZ 

0.769 

3.004 

c 3 h 

aug-cc-pVDZ 

cc-pVDZ 

- 

3.512 


aug-cc-pVTZ 

cc-pVTZ 

- 

3.551 

c 4 h 

aug-cc-pVDZ 

cc-pVDZ 

0.870 

4.404 

C«jH 

aug-cc-pVDZ 

cc-pVDZ 

- 

4.881 

c 6 h 

aug-cc-pVDZ 

cc-pVDZ 

0.979 

5.536 

c 7 h 

aug-cc-pVDZ 

cc-pVDZ 

- 

5.945 


be due to the vibration effects incorporated in the 
latter value. It is noteworthy that B c is very sensitive 
to the C-C bond length: a change of only 0.001 A 
shifts B e by 60 MHz. In fact, including core and 
core-valence correlation increases the error by more 
than 200 MHz. 

The RCCSEKT) dipole moments /x e of C„H radi- 
cals are given in Table 2. It is evident from C 2 H that 
Me is significantly changed by adding diffuse func- 
tions to C, but not to H. Therefore, this pairing of the 
aug-cc-pVXZ set for C with the standard cc-pVXZ 
set for H was used for the larger species as well. 
Furthermore, the results for C 2 H and C 3 H indicate 
that there is little to be gained by improving from 
double to triple zeta quality sets. The impact is 
negligible for the 2 S + state of C 2 H and only about 
0.05 D for the 2 II states of C 2 H and C 3 H. The 
correlated RCCSEKT) values are expected to be 
somewhat more reliable than the RHF or UHF val- 
ues reported in previous work [1,3]; unfortunately, 
there are no experimental values for comparison. It is 
of interest to note that although the dipole moments 
of both states increase with chain length, the magni- 
tude doubles for the 2 II state from C 2 H to C 7 H. Of 
note: although a large dipole moment improves the 
likelihood of detecting a species, it also increases the 
long range attraction for other molecules and may 
lead to rapid destruction of the longer C„H radicals 
that have not been detected. 
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4. The 2 %+- 2 U state separations of C 2 H, C 4 H, 
and C 6 H 

The determination of accurate excitation energies 
requires a balanced treatment of the two states with 
respect to both methodology and the one-electron 
basis sets employed. In order to judge the expected 
accuracy of RCCSD(T)/cc-pVXZ calculations for 
describing the "2 + - 2 n equilibrium state separation 
T t of C n H species, we will begin by benchmarking 
the behavior of C 2 H. Although there have been 
several attempts to measure this quantity, a high 
accuracy value has not been determined. Values have 
ranged from 1600 to 3800 cm' 1 [25]. The best 
theoretical work on C 2 H has been reported by Krae- 
mer et al. [26] ( T c = 3650 cm' 1 ), Koures and Hard- 
ing [27] (T e = 3550 cm -1 ), Thummel et al. [28] 
(T c = 3300 cm -1 ) and Peric et al. [29], which should 
be consulted for citations to other previous ab initio 
work. 

Table 3 lists total energies and T c for CMRCI, 
CMRCI + Q, and RCCSD(T) calculations on C 2 H 
with the cc-pVDZ, cc-pVTZ, and cc-pVQZ basis 
sets and the extrapolation to the CBS limit. Given 
that it is not unusual to find CMRCI and CMRCI 4- Q 
calculations bracketing experiment [30], the interme- 
diate location of the RCCSD(T) result lends credibil- 
ity to using the method for C 4 H and C 6 H. 


Table 3 

Total energies E e and excitation energy T c of the 2 X + and 2 n 
states of C 2 H at the indicated levels of theory 


Method 

Basis set 

E c (£ h ) 
2 £ + 

2 n 

T e (cm 1 ) 

CMRCI 

cc-pVDZ 

- 76.39842 

-76.38398 

3169.2 


cc-pVTZ 

-76.46115 

-76.44492 

3562.1 


cc-pVQZ 

-76.47958 

-76.46280 

3682.8 


est. CBS limit 

-76.48987 

-76.47277 

3753 

CMRCI + 0 cc-pVDZ 

-76.40374 

-76.39059 

2886.1 


cc-pVTZ 

-76.46961 

-76.45517 

3169.2 


cc-pVQZ 

-76.48892 

-76.47405 

3263.6 


est. CBS limit 

- 76.49969 

-76.48458 

3316 

RCCSIXT) 

cc-pVDZ 

- 76.40063 

-76.38738 

2908.0 


cc-pVTZ 

-76.46779 

-76.45301 

3243.8 


cc-pVQZ 

-76.48763 

-76.47236 

3351.4 


est. CBS limit 

-76.49872 

-76.48316 

3415 


cc-pCVTZ 

-76.569561 

-76.554047 3404.9 


Table 4 

Total energies E c and excitation energy T t of the "3C (even n) 
and 2 n states of C„H at the RCCSD(T) level; a positive value of 
T c indicates that the state is more stable 


Species 

Basis set 

E'(E h ) 


(cm 1 ) 


2 n 

C,H 

cc-pVDZ 

- 

-114.35711 



cc-pVTZ 

- 

-114.45638 



cc-pVQZ 

~ 

- 114.48575 



est. CBS limit 

- 

- 114.50216 


c 4 h 

cc-pVDZ 

- 152.33985 

-152.34105 

-263 


cc-pVTZ 

- 152.47463 

- 152.47487 

-53 


cc-pVQZ a 

- 152.50249 

- 152.50236 

29 


est. CBS limit 

- 152.51664 

- 152.51631 

72 

c 5 h 

cc-pVDZ 

- 

- 190.30893 



cc-pVTZ 

- 

- 190.47665 


c 6 h 

cc-pVDZ 

- 228.28273 

- 228.28764 

-1078 


cc-pVTZ a 

-228.48539 

-228.48946 

-895 

c 7 h 

cc-pVDZ 

- 

- 266.25606 



11 Single point calculation was performed using estimated bond 
lengths (see text). 


The computed total energies for C 3 H through 
C 7 H as well as the T e values for C 4 H and C 6 H are 
in Table 4. For C 4 H, the most recent correlated 
calculations for T c are the combined experimental 
and theoretical study reported by McCarthy [7] and 
the work of Natterer and Koch [8]. Since the former 
study is quite similar to the present one, this discus- 
sion will focus on the latter work. Natterer and Koch 
found the 2 £ + state to be more stable at the MCSCF 
and MRCI levels, but the 2 n state was more stable 
at the MRCI + Q, ACPF, and RCCSD(T) levels of 
theory. For the ACPF and MRCI + Q methods, the 
2 n state was more stable than the 2 £ + state by 100 
to 240 cm' 1 , depending upon the basis set used. On 
the other hand, T t jumped to — 1 1 37 cm ~ 1 or more 
for the RCCSEXT) method. The present calculations 
are at odds with these results. With the smallest 
cc-pVDZ sets, the 2 II state is more stable, but only 
by 263 cm -1 . Improving the basis sets lowers the 
state preferentially, and it becomes more stable 
at the cc-pVQZ level, though only by 29 cm' 1 . 
Extrapolating to the CBS limit stabilizes the 2 X + 
state to 72 cm' 1 relative to the 2 31 state. The 
benchmark calculations on C 2 H suggest that, if any- 
thing, RCCSD(T) may be underestimating T e . Unfor- 
tunately, Natterer and Koch did not report calcula- 
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tions on the state separation in C 2 H, which would 
have provided a test on the approach they adopted 
for C 4 H. Aside from using different basis sets, they 
also ran the correlated calculations at their MCSCF 
optimized geometries. Their bond lengths are roughly 
comparable to the present RCCSD(T)/cc-pVTZ val- 
ues, so this cannot account for the 1000 cm -1 differ- 
ence in T c at either the cc-pVTZ level. 

Natterer and Koch were forced by convergence 
problems in their RHF and UHF calculations to 
adopt the quasi-RHF (QRHF) method to generate 
zeroth-order wave functions. When MOLPRO was 
used to run an RCCSEKT) calculation based on an 
RHF wave function with the [4321/321] ANO basis 
set (yielding respective total energies for the 2 FI and 
2 S + states of -152.47037 and -152.46864 £ h ), 
the 2 I1 state was found to be only 378 cm 1 more 
stable than the 2 S + state, rather than the 1137 cm 5 
reported in Ref. [8]. The new value is more consis- 
tent with the present results and implies that the 
QRHF wave function introduces a large amount of 
error. Also, the MRCI and ACPF results noted above 
are closer to the RHF-based value of T t . 

While the cc-pVTZ and [4321/321] ANO sets 
yield similar total energies, their respective T e values 
differ by more than 300 cm -1 . In order to test the 
accuracy of the ANO sets for describing state separa- 
tions, the [4321 /321] ANO contraction was used to 
determine T e for C 2 H (again at the cc-pVTZ geome- 
try). A value of 2785 cm ' 1 was obtained (with total 
energies for the 2 2 + and 2 n states of -76,46284 
and —76.45015 E h \ which is about 460 cm -1 
below the cc-pVTZ result and in much poorer agree- 
ment with typical theoretical values. The trend for 
both C 2 H and C 4 H indicates that the [4321/321] 
ANO set performs less accurately than the cc-pVDZ 
set. The ANO sets are also much more computation- 
ally demanding. The present conclusion that the 
ground state of C 4 H is a 2 X* state is reinforced, 
although clearly the 2 11 state must lie remarkably 
close. 

For C 6 H, previous correlated efforts include the 
Cl work of Pauzat and Ellinger [2], where the 2 n 
state was found to be more stable than the 2 X + state 
by 6385 and 2379 cm -1 in first-order (FOCI) and 
limited second-order (SOC1) calculations, respec- 
tively. Natterer and Koch [8] reported a state separa- 
tion of 2175 cm -1 at the QRHF + CCSD(T)/ 


[4321/321] level. The best present value, using the 
cc-pVTZ sets and estimated geometries, places the 
2 n state 895 cm -1 below the 2 S + state. 

Core and core-valence correlation impact excita- 
tion energies as well as bond lengths. Table 3 shows 
that the state of C 2 H is stabilized by an addi- 
tional 160 cm 1 with respect to the valence-electron 
triple zeta result when the cc-pCVTZ sets were used 
in all-electron calculations. If the trend observed in 
CH 2 [12] carries over to C 2 H, this correction may be 
too small by about 50% with respect to the CVXZ 
limit. The shift in C 2 H suggests that the state separa- 
tion in C 4 H may be somewhat larger than the <100 
cm -1 valence CBS limit, perhaps reaching the 200- 
300 cm' 1 range. 

5. Electron affinities of C 2 H and C 3 H 

The electron affinities of C 2 H and C 3 H were 
computed at the RCCSEKT) level by optimizing 
C 2 H" and C 3 H' with the aug-cc-pVTZ set on C 
and the cc-pVTZ set on H, which required reoptimiz- 
ing the neutrals with the same sets. Geometries, 
energies and EAs are presented in Table 5. This 
strategy was sufficient to reproduce the EA of C 2 H 
to within the experimental uncertainty (2.969 ± 0.006 
eV, [31]), although the inclusion of a partial zero- 

Table 5 

Structures (r c , A; angle, deg) and total energies (£ e , £ 0 , E h ) of 
C 2 H, C 3 H and their anions and adiabatic electron affinities (EA, 
eV) at the RCCSD(T)/{C: aug-cc-pVTZ; H: cc-pVTZ} level; 
experimental values are C 2 H (2.969 ± 0.006) and C 3 H (1.858 ± 
0.023) 


Species 

Quantity 

Neutral 

Anion 

C 2 H 

r e (C,-C,) 

1.2152 

1.2549 


r e (C 2 -H) 

1.0658 

1.0726 



-76.47122 

-76.58015 


(EA) C 

2.964. 



E o 

-76.45727 

-76.56592 


(EA)„ 

2.957 


c 3 h 

r e (C,-C 2 ) 

1.3468 

1.2897 


r e (C 2 -C 3 ) 

1.2484 

1.3605 


r e (C,-H) 

1.0654 

1.1142 


zKC.-Cj-Cj) 

180 

174.73 


^(C 2 -C 3 -H) 

180 

- 108.25 


E c 

-114.46241 

-114.52949 


(EA) C 

1.825 
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point correction (omitting the frequencies associated 
with the Renner-Teller split bending modes) in- 
creases the error slightly, to 0.012 eV. Values of 
3.15 and 2.96 eV, computed by Lima and Canuto 
[32] and Natterer and Koch [8], respectively, have 
previously been reported. Both C-C and C-H bonds 
length increase slightly in the anion, a state, but 
the effect is larger ( ~ 0.04 A) for the C-C bond. 

While C 2 H is linear, C 3 H is not. The ground 
state is of f A' symmetry, and there is a low-lying, 
bound excited state of 3 A" symmetry. For the *A 
state, the C-C-C angle remains quasilinear ( = 175°), 
but the C-C-H angle has decreased to 141.2°. The 
C-C-C-H dihedral angle is 180°, making the termi- 
nal C and H trans with respect to the central C pair. 
The two C-C bond lengths change significantly with 
respect to the neutral. The computed equlibrium EA 
is 1.825 eV, which falls just slightly below the 
experimental range of (1.858 ± 0.023) of Oakes and 
Ellison [33]. Our result appears to be the first theoret- 
ical prediction of this quantity. The 3 A" state (E e = 
-114.52730 E h ) lies only 481 cm 1 above the 
ground state with the triple zeta sets, but its geome- 
try is somewhat different: /* e (C 1 -C 2 ) = 1.309 A, 
r e (C,-C 3 )= 1.312 A, r c (C 3 -H)= 1.076 A, Z(C- 
C-C)= 174.66° and Z(C-C-H) = 141.20°. Dua et 
al. [9] reported structures optimized at the HF level 
for both states that are similar to the RCCSEXT) 
ones, but they found the 3 A" state to be slightly 
nonplanar. They also reported the triplet state to be 
lower in energy, by 7050 cm 1 at the HF/6-311 + 
+ G* * level and by 1236 cm 1 at the CID/6-311 
+ + G * V/HF/6-311 + -F G * * level. 

6. Conclusions 

The RCCSIXT) method and correlation consistent 
basis sets were used to determine the structures and 
other properties of C„H radicals for n = 2-7. In 
general, predicted valence-correlation values of the 
rotational constants B e were within 100 MHz of 
experimental or astronomical B 0 values; adding a 
correction for core and core-valence contraction im- 
proves the agreement substantially. The predicted 
S + - 2 n state separation of C 2 H at the RCCSD(T) 
level of theory is in good agreement with values 
determined using multireference methods and with 


other theoretical results. In the estimated complete 
basis set limit, the ground states of C 2 H and C 4 H 
were found to be of symmetry and lie 3415 and 
72 cm 1 below their respective 2 IT states. The Yl 
ground state of C 6 H lies 895 cm 1 below the + 
state (RCCSD(T)/cc-pVTZ). 

The electron affinities of C 2 H and C 3 H were 
determined to be 2.964 and 1.825 eV, respectively, 
which agree well with the experimental values of 
2.969 + 0.006 and 1.858 + 0.023 eV. Although 
C 2 H is linear, C 3 H is not: its ground state is a 
planar *A' state. A 3 A" state lies only 481 cm 1 
higher in energy. 
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